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XVII. On the Secular Variation of the Moon's Mean Motion. 
By J. C. Adams, Esq.,M.J., F.R.S. 8fc. 

Received and Read June 16, 1853. 

1. IN treating a great problem of approximation, such as that presented to us by 
the investigation of the moon's motion, experience shows that nothing is more easy 
than to neglect, as insignificant, considerations which ultimately prove to be of the 
greatest importance. One instance of this occurs with reference to the secular 
acceleration of the moon's mean motion. Although this acceleration, and the dimi- 
nution of the eccentricity of the earth's orbit, on which it depends, had been made 
known by observation as separate facts, yet many of the first geometers altogether 
failed to trace any connexion between them, and it was only after making repeated 
attempts to explain the phenomenon by other means, that Laplace himself succeeded 
in referring it to its true cause. 

2. The accurate determination of the amount of the acceleration is a matter of very 
great importance. The effect of an error in any of the periodic inequalities upon 
the moon's place, is always confined within certain limits, and takes place alternately 
in opposite directions within very moderate intervals of time, whereas the effect of an 
error in the acceleration goes on increasing for an almost indefinite period, so that 
the calculation of the moon's place for a very distant epoch, such as that of the 
eclipse of Thales, may be seriously vitiated by it. 

In the Mecanique Celeste, the approximation to the value of the acceleration is 
confined to the principal term, but in the theories of Damoiseau and Plana the 
developments are carried to an immense extent, particularly in the latter, where the 
multiplier of the change in the square of the eccentricity of the earth's orbit, which 
occurs in the expression of the secular acceleration, is developed to terms of the 
seventh order. 

As these theories agree in principle, and only differ slightly in the numerical value 
which they assign to the acceleration, and as they passed under the examination of 
Laplace, with especial reference to this subject, it might be supposed that at most 
only some small numerical corrections would be required in order to obtain a very 
exact determination of the amount of this acceleration. 

It has therefore not been without some surprise, that I have lately found that 
Laplace's explanation of the phenomenon in question is essentially incomplete, and 
that the numerical results of Damoiseau's and Plana's theories, with reference to it, 
consequently require to be very sensibly altered. 
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3. Laplace's explanation may be briefly stated as follows. He shows that the mean 
central disturbing force of the sun, by which the moon's gravity towards the earth is 
diminished, depends not only on the sun's mean distance, but also on the eccentricity 
of the earth's orbit. Now this eccentricity is at present, and for many ages has been, 
diminishing, while the mean distance remains unaltered. In consequence of this 
the mean disturbing force is also diminishing, and therefore the moon's gravity 
towards the earth at a given distance is, on the whole, increasing. Also, the area 
described in a given time by the moon about the earth is not affected by this alteration 
of the central force ; whence it readily follows that the moon's mean distance from the 
earth will be diminished in the same ratio as the force at a given distance is increased, 
and that the mean angular motion will be increased in double the same ratio. 

4. This is the main principle of Laplace's analytical method, in which he is followed 
by Damoiseau and Plana ; but it will be observed, that this reasoning supposes that 
the area described by the moon in a given time is not permanently altered, or in 
other words, that the tangential disturbing force produces no permanent effect. On 
examination, however, it will be found that this is not strictly true, and I ivill endeavour 
briefly to point out the manner in which the inequalities of the moon's motion are 
modified by a gradual change of the central disturbing force, so as to give rise to 
such an alteration of the areal velocity. 

As an example, I will take the Variation, the most direct effect of the disturbing 
force. 

In the ordinary theory, the orbit of the moon as affected by this inequality only, 
would be symmetrical with respect to the line of conjunction with the sun, and the 
areal velocity generated while the moon was moving from quadrature to syzygy, 
would be exactly destroyed while it was moving from syzygy to quadrature, so that 
no permanent alteration of areal velocity would be produced. 

In reality, however, the magnitude of the disturbing force by which this inequality 
is caused, depends in some degree on the eccentricity of the earth's orbit, and as this 
is continually diminishing, the central disturbing forces at equal angular distances 
on opposite sides of conjunction will not be exactly equal. Hence the orbit will no 
longer be symmetrically situated with respect to the line of conjunction. Now the 
change of areal velocity produced by the tangential force at any point, depends partly on 
the value of the radius vector at that point, and consequently the effects of the tangential 
force before and after conjunction will no longer exactly balance each other. 

The other inequalities of the moon's motion will be similarly modified, especially 
those which depend, more directly, on the eccentricity of the earth's orbit, so that 
each of them gives rise to an uncompensated change of the areal velocity. 

Since the distortion in the form of the orbit just pointed out is due to the alter- 
ation of the disturbing force consequent upon a change in the eccentricity of the 
earth's orbit, and it is by virtue of this distortion that the tangential force produces 
a permanent change in the rate of description of areas, it follows that this alteration 
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of the areal velocity will be of the order of the square of the disturbing force multi- 
plied by the rate of change of the earth's eccentricity. 

It is evident that the amount of the acceleration of the moon's mean motion will 
be directly affected by this alteration of areal velocity. 

5. Having thus brieflyindicated the way in which the effect now treated of originates, 
I will proceed with the analytical investigation of its amount. 

In the present communication, however, I shall confine my attention to the prin- 
cipal term of the change thus produced in the acceleration of the moon's motion, 
deferring to another, though I hope not a distant, opportunity, the fuller development 
of this subject, as well as the consideration of the secular variations of the other 
elements of the moon's orbit arising from the same cause. 

In what follows, the notation, except when otherwise explained, is the same as that 
of Damoiseau's " Th^orie de la Lune." 

6. If we suppose the moon to move in the plane of the ecliptic, and also neglect the 
terms depending on the sun's parallax, the differential equations of the moon's motion 
become 

° = dF+«-T.+ 2tf*+2 IV COS (2,-2,') 

3 m'u' 3 du . , n „„ 3m'/ ., d*u\Cv! 3 dv . .„ „ 

-2AV S - 8,n ( 2 '- 2 ')-ir(« + ^;J- 7 - «n(2i-20 

dt 1,3 m> Cu'Hv . . nl . ,27 m' 2 rCu' 3 d V . , n n . ' 

In the solution usually given of these equations, u is expressed by means of a 
constant part and a series involving cosines of angles composed of multiples of 
2v— 2mv, cv—m, and c'mv — m'; also t is expressed by means of a part proportional 
to v and a series involving sines of the same angles ; the coefficients of the periodic 
terms being functions of m, e and e'. Now if e' be a constant quantity, this is the 
true form of the solution, but if e' be variable, it is impossible to satisfy the differen- 
tial equations without adding to the expression for u a series of small supplementary 
terms depending on the sines of the angles whose cosines are already involved in it, 
and to that for t, similar terms depending on the cosines of the same angles, the 

coefficients of these new terms involving -r. as a factor. 

The quantity l-^r sin (2v— 2/), which occurs in the above equations, is proportional 

to the variable part of the square of the areal velocity, and consists, in the ordinary 
theory, of a series of periodic terms involving cosines of the angles above mentioned. 
In consequence, however, of the existence of the new terms just described, there will 
be added to it a series of small terms involving sines of the same angles, together 

with a non-periodic part of the form jHe'rfe' or ^ He' 2 . The introduction of this term 
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will evidently change the relation between the non-periodic part of -j and e' 2 , upon 

which the secular acceleration depends. 

7. We must commence by finding the new terms to be added to the ordinary 
expression for u. 

For the sake of simplification we will neglect the eccentricity of the moon's orbit. 

Let - denote the non-periodic part of u, and — \-lu the complete value. 

Then by substitution in the equation for u, making use of Damoiseau's develop- 
ments of the undisturbed values of the several functions of u, u', and c— c' which occur 

in it, putting A 2 =a /5 and writing, for convenience, mv instead of 1 mdv-\-\, and c'mv 
instead of c' I mdv-\-X— w' (as in Plana, vol. i. p. 322), we obtain 

n d \a) 1 1 ,d*$u,. 

+2 TV+arJ+H a ° B +2 V e c°S cmv — 2 7-0 +3e cos c'mi>}a5u 

d(-) 
B —„ a \aJ . ,_ „ . , 3 m 2 / 5 ,A .„ 

— 2 mr J -81 " (2v — 2wij') + 2 — ( 1— ^e J cos (2)'— 2mc) 

. 21 jn 2 , ._ , . 3 to 2 . , . 

+^7 — e cos (2c — 2mv — cm») — j — e cos (2c — 2mv-\-cmv) 

3m' 



— l^n( 1 — 2 e ' 2 J s i n (2v — 2?»v)+2e f sin (2c— 2mv— c'mv) 

— ^e' sin (2c— 2mc-j-cW) 
— j ( 1— 2 e y cos {^—2mv)-\-^ cos (2c— 2mc— c'mv) 



— ~e'cos (2c — 2mv-\-dmv) \atu 

1 * J 

3 m 2 f / 5 \ 7 

— 2 — ( 1 — 2 e ' 2 J sm (2"— 2mc) +7: e' sin (2c— 2mc— c'mc) 

— ^e' sin (2c— 2mc+c'mc) f-^r- 

+ 12 — \dv\ i 1— ^e' 2 ) sin (2c— 2mv)-\--d sin (2c— 2mc— c'mv) 

— g^'sin (2c— 2mc-f-c'mc) >«5w 

-^l^^+^jWf 1 -^' 2 ) sin (2c-2mc)+|e , sin(2c-2»ic- C 'r 

— se' sin (2c— 2m»-rc'm») •. 



■ )»» 



OF THE MOON'S MEAN MOTION. 401 

8. Also, assume 

(5 \ dde* 

1— jje'M cos (2c— 2»ii')+a 30 ^-sin (2v— 2mv) 

—-^rrire cos cmv+a, 6 -^ sin cmv 

-J-h»wV cos (2v—2mv—c'mv)-\-a 33 -j.s\a (2v—2mv — c'mv) 

1 de' 

— gwVcos (2»— 2mi>+c'mv)+a 3i ^j t sin (2v—2mv-\-c'mv), 

where the coefficients of the terras involving cosines are those given by the ordinary 
theory, and a 30 , a l6 , a 33 , and a 3i are numerical quantities to be determined. 

9. In developing the terms of the above equation, by the substitution of this value 

dd dd 

of alu, the quantity ^ may be considered constant, and -r must be expressed in terms 

of it. 

,_, de! ndt dd 

lhUS Tv—lh^dt 

de' f ] 1 77 

~~dt\ l~T m2cos (2 v —2 mv ) — g-mV cos (2v— 2mv— c'mv) 

+-g-mV cos (2f— 2mv-\-c'mv) I. 

Also, integrating by parts, and putting 2 instead of 2— 2m, 2— '6m, and 2— m in 
the divisors introduced by integration, since we only want to find the terms of the 

dd 
lowest order which are multiplied by ^-, we obtain 

3m 2 r f/ 5 \ 7 

\dv\ ( 1— 2 e ' 2 ) sni (2c— 2mc)+2 e ' sm (2v— 2mc— c'mc) 

— ^e' sin (2?— -2»w+c'»»') f 

3 »i 2 / 5 \ 21 ?n 2 

= ~ — ( 1 — -e' 2 ) cos (2v—2mv)-\--r — d cos (2v— 2mv— dmv) 

" a i\ 2 J 4 «; 

3 m 2 
— t — e' cos {2v— 2mv-\-dmv) 

, 15 m*C, e'dd ndt , n „ . 21 m 2 r, rfe 7 ra<ft , n n , . 

+¥ #5# -*T C0S (2<'-2mv)- T -Jcfi,^ ^-cos (2»-2itw-c'm») 

, 3 m 2 f , <fe' »wft /««.»» 
+4 Ty^t ~dJ cos (2"-2mc+c'mv). 

And aW=3mV sin cW[— e* sin <,"'mv] 

= - 2 -mY, 

retaining only the term which will be required. 

3 62 
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10. When the proper substitutions are made, the terms involving cosines destroy 
each other, as in the usual theory, and by equating to zero the terms involving the 

sines, we obtain 

15 
20m 2 — 3a m -\--jm 2 =0, 

95 2 95 2 

or 3a S0 =-pn 2 :. a 30 =— m 



3m 3 +a l6 =0 :. a l6 = — 3m 3 

21 

T 



21 
— 14m 2 — 3a 33 — — m 2 =0, 



or aa„= — T-m° :. «„„= — tttw? 



or 



"^33 — 


133 

:-^-m 2 


.•. a. 


J3 


2m 2 - 


3 

•3a 3 4+^p 2 


=0, 




3«34 = 


19 2 . 


«34 = 


19 2 

: 24 m 



133 

24 



11. In order to obtain the relation between a and a p we must substitute the value 
just found for alu, in the same equation, and equate to zero the non-periodic part, 
observing that the terms 

12^ f«frj(l — |e' 2 ) sin (2»— 2»w) + |e' sin (2k— 2rm-c'mv) 

— xe'sin (2v—2mi/-\-c'mi>) \ahu 



give 



12m 2 C j f 95 /*' 931 /*' 19 *<?<&) 
285 m 4 f 7 eW 
285 m 4 ,„ . . ... 

= — g- — e 2 as their non- periodic part. 
Also the terms 

15 m 2 f , eW ndt 21 m z f , ^ »<# /no r \ 

+4 ^ J*^ A" cos ( 2 "~ 2^"+cW) 
of art. 9. similarly give 

15 m 2 C,f 11 2 *\ 21m 2 f , / 77 /&'\ ,3 m 2 f ,/U /^A 

_ 165 ro 4 1617 m 4 , 2 33 m 4 „ 

= -l2"^ e + 12§V +128 ^ e near, y 

495 m 4 

=-gj — e' 2 as their non-periodic part. 
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i 2. Hence we obtain 



0=1-1+1™ (i+f^H^+f ^+f ^ 

a «, 2 «, V 2 / 4 «, ' 64 «, '8 a, 



9 »t' 



4 « . \ / 



441 w 4 ,„ 9 w 4 



4 a t 

3 m 4 
'2 « 



-e' 2 -, 



16 a t " 16 a y 



+;:-(.-^>+^>+i">- 2 ?=v 



8 «, 



8 a, 



8 a. 



4 // v 1 — i}e )— Ifi „ e — |« „ ^ J 



16 a, 16 a, 



or 



llf 1 _ ^ _ ^l'i 5 ? 1 



Now m 2 = ( t ■ „\3 in Piana's notation, or (substituting the value of p given in Plana, 
vol. ii. p. 855). 

m 2 =m 2 (1 —pn i —-m i e! 2 J nearly, 



.1 1/ 1 2 13 4 3 2 , 2 .3201 4 , 2 

.. -=H1- 7;m 2 +-T-m 4 — T»reN — — m 4 e 

a a.\ 2 ' 4 4 ' 



64 



and 



a 2 =a, 2 { 



, |_ • 23 4 , 3 2 12 3129 • 



V s ] . 



]3. Again, by substitution in the equation for -v-, we obtain 
|=-^ 1 -2a&w+|m 4 (l -5e' 2 ) + | 7 mV 2 +^mV 2 +gmV 2 

' 1 jM(>- 



! 3 "I 8 

+ 2»» 2 - 



5 \ 7 

ge' 2 J sin (2v— 2mv)-\- „-e' sin (2v— 2»jv— c'mv) 



— ~e' sin (2i>— 2wi'+cW) 



— 3m 2 ~ alu \dvU\—^\ sin (2v— 2»w) + |e' sin {2v—2mv—chnv) 

—d sin (2f— 2mv+c!mv) 
— 6rf^ j ^[(l— 2 e ' 2 ) s i n ( 2 "— 2mv) + |e'sin (2k— 2ms<— dmv) 

— ge'sin (2v— 2rm-\-dmv) \a$u 

27 /a\ 2 ( C f f 5 \ 7 

+-g-m 4 ( -J j i e?e ( 1— xe' 2 ] sin (2v— 2mv)-\-^d sin (2f— 2mv— c'mv) 

—-xd sin (2f— 2/mf+c'mj») \ . 

14. Develope this equation as before, retaining m 4 only when it occurs in the 
non-periodic part, and we have 
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S=-^ 1 -2alu+\m*+\m\ 1 -5 e ' 2 )+|m*(l -Se' 2 )-^' 2 

— 7»ra 2 ( 1 — 2 e ' 2 ) cos (2c— 2mc) — -g-mV cos {2v—2mv—c'mv) 



3 

+g»?y cos (2c— 2mc+cW) 

15 e'rfe' 21 de 1 

■-^m 2 -jTs\n (2v—2mv) -^YQm 2 ^ sin {2v—2mv—dmv) 



— YgW 2 ^sin (2c— 2mv+dmv) L 



d* a 2 f, ,171 . , 2391 , „, 



l / 5 \ 425 ede 

-mHl—^e' 2 ) cos (2c— 2mc)— -gjw^sin (2c— 2mc) 



de' 
+37wV cos c'mv-\-6m a —r t sin cW 

7„_ 2 _. ,_ „ . v , 595 „de' 



—~m 2 e' cos (2c— 2mc— c'mc) +^-»i 2 ^ sin (2c— 2mc— cmc) 
+-mV cos (2c— 2mc+c'wjc)— jgWi 2 ^ sin (2c— 2mc+c'mc) L 

15. Substitute the value before found for a 2 in terms of a, 2 ; 

-- 7m 2 [l-^' 3 j COS (2c— 2?«c) — 04 m2 ^ sm (2c— 2mc) 

de' 
+3mV cos c'mc+6m 3 -^ sin cW 

77 595 de' 

— -rrmV cos (2c— 2»ic — c'mc) + ^rm 2 -^ sin (2c— 2mc— cmc) 

+-j-mV cos (2c— 2mc+cW)— ^g-m 2 ^- sin (2c— 2mc+c'mc) I. 

16. Now, put i=af|l +m 2 -^+^' 2 -^-Wj, 
multiply by w, and integrate ; 

J 11 / 5 \ 295 e*de' 

ndt=v—-^m 2 ( 1 — ^e' 2 J sin (2c— ^m»)-\--^m 2 ^ cos (2c— 2mc) 



de' 



mc 



+3zrae' sin c'mc-f3—T- cos c' 

77 413 de' 

— jgmV sin (2c— 2mc— cW)— ^g m 2 ^cos (2c— 2mc— cW) 

11 59 de 1 

+Yg»iV sin (2c— 2mc+c'wc)+^gm 2 ^ cos (2c— 2mv+cmv). 
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17. In the expression for -just found, a t is absolutely constant, but e' is variable, 

consequently n will vary, and therefore m likewise, which is connected with it by 

n' 
the equation m=-- 

Taking the variation of the equation for n, and observing that — = — — , we have 
o=^ (1 -„.) + (| K ._5|fV>( e »), 

Therefore, if N be the initial value of n, and E' the corresponding value of e', 

n:= N-(fm 2 -^V)ra(e' 2 -E' 2 ), 

and (Lft=N*+s-(|m 2 -^m 4 )(V 2 -E'><ft. 

Hence the expression for the true longitude in terms of the mean, contains the secular 
equation 

18. According to Plana, the corresponding terms in the expression for the secular 
equation are 

-(f™ 2 -^™ 4 )J(e' 2 -E><ft. 

Hence we see that the terms now taken into consideration have the effect of making 
the second term of the secular equation more than three times as great as it would 
otherwise be. Of course, the succeeding terms will also be materially changed. 
The principal term of the correction to be applied to Plana's value of the secular 



acceleration is therefore 

5355 



m 4 f(e' 2 -E'>^. 



128 
Now f(e' 2 -E><ft= - 1270"(~) 2 nearly, 

where t is expressed in years ; therefore the numerical value of this term is 

- 1 "- 66 (4) 2 ' 

This result will serve to give an idea of the numerical importance of the new terms 
to be added to the received value of the secular acceleration, and probably will not 
differ widely from the complete correction; though in order to obtain a value suffi- 
ciently accurate to be definitively used in the calculation of ancient eclipses, the 
approximation must be carried considerably further. 

The new periodic terms added to the moon's longitude are perfectly insignificant, 
the coefficient of that involving cos c'mv, which is by far the largest of them, only 
amounting to 0"'003. 
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19. Transforming the expressions found above, so as to obtain the moon's longitude 
and radius vector in terms of the time, and writing for convenience nt instead of 

Jndt-\-s, mnt instead of mnt-\-s\ and c'mnt instead of c'mnt-\-i'— w', we have 

11 / 5 \ 74 e'del 

v=nt+— m 2 ll— ge'M sin (2 — 2m)nt— y?» 2 ^ cos (2— 2m)nt 

— dme sin cmnt — 3—r. cos cmnt 

+ jg?wV sin (2— 2m — c'm)nt-\- ~7Q m2 ~jt cos (2—2m—c'm)nt 

11 257 a 

— jg/raVsin (2— 2m-\-c'm)nt— "lg m2 ^ cos (2— 2m+c'm)n£ 

o , 11 . 201 , ,„ 

- = au — 1 — T.~m — -rx-mre 2 
r 8 16 

(5 \ 203 e'de* 

l ~2 e 'j C0S ( 2 ~ 2m ) nt ~^^2 m2 ^di Sin ( 2 — 2m ) nf 

3 <fe' 

-^mV cos c'mwi — 3m 3 — r. sin c'wwtf 

7 61 de* 

+ o»tVcos (2— 2m— c'm)nt— ^m^sin (2—2m—c'm)nt 

1 91 rfe* 

— 2»J.Vcos (2— 2m+c'm)n^+^m 2 -T>sin (2 — 2m+c'm)wf. 

20. The existence of the new terms in the expressions for the moon's coordinates 
occurred to me some time since, when I was engaged in thinking over a new method 
of treating the lunar theory, though I did not then perceive their important bearing 
on the value of the secular equation. 

My attention was first directed to this latter subject while endeavouring to supply an 
omission in the theory of the moon given by Pontecoulant in his " Th6orie Analy- 
tique." In this valuable work, the author, following the example originally set by 
Sir J. Lubbock in his Tracts on the Lunar Theory, obtains directly the expressions 
for the moon's coordinates in terms of the time, which are found in Plana's theory 
by means of the reversion of series. With respect to the secular acceleration of the 
mean motion, however, Pontecoulant unfortunately adopts Plana's result without 
examination. On performing the calculation requisite to complete this part of the 
theory, I was surprised to find that the second term of the expression for the secular 
acceleration thus obtained, not only differed totally in magnitude from the correspond- 
ing term given by Plana, but was even of a contrary sign. My previous researches, 
however, immediately led me to suspect what was the origin of this discordance, and 
when both processes were corrected by taking into account the new terms whose 
existence I had already recognized, I had the satisfaction of finding a perfect agree- 
ment between the results. 

June \6th, 1853. 



